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Fully frustrated Josephson Junction arrays (FF-JJA's) exhibit a subtle compound phase tran- 
sition in which an Ising transition associated with discrete broken translational symmetry and a 
Berezinskii-Kosterlitz-Thouless (BKT) transition associated with quasi-long-range phase coherence 
occur nearly simultaneously. In this Letter we discuss a cold atom realization of the FF-JJA sys- 
tem. We demonstrate that both orders can be studied by standard momentum-distribution-function 
measurements and present numerical results, based on a successful self-consistent spin- wave approx- 
imation, that illustrate the expected behavior of observables. 

PACS numbers: 03.75.Lm 
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The preparation of cold atomic gases trapped in an 
optical lattice has opened up attractive new possibil- 
ities for the experimental study of strongly correlated 
many-particle systems and has inspired much the- 
oretical activity (see e.g. Ref. 2 for a review). In 
particular, the experimental observation by Greiner et 
al. of a Superfluid-Mott Insulator (SI) transition in a 
three-dimensional (3D) optical lattice explicitly demon- 
strated the possibility of realizing strongly-correlated 
cold bosons. The SI transition in an optical lattice was 
predicted in Ref.jS and can be described by the Bose- 
Hubbard model which has also been employed to 
model 2D granular superconductors 5] and JJA's 0. 
This success has motivated many new proposals for 
cold-atom simulations of strongly correlated boson phe- 
nomena. 

In this Letter we propose that cold atoms be used to 
study the incompletely understood phase transitions that 
occur in FF-JJA's 0,0. The boson Hubbard model for 
JJA's accounts for Cooper pair hopping between small 
superconducting particles and for Coulomb interactions 
which can be dominantly intra-particle. For supercon- 
ducting particles the model applies when the thermal en- 
ergy ksT is much smaller than the bulk energy gap, i.e. 
when the underlying fermionic character of electrons is 
suppressed. Cold atoms in optical lattice potentialspro- 
vide, in some senses at least, a closer realization 0, of 
the boson Hubbard model because other degrees of free- 
dom are more completely suppressed and because the 
interactions are more dominantly on-site. Frustration 
can be introduced into JJA's by introducing an external 
magnetic field to change the energetically preferred phase 
relationship between boson amplitudes on neighboring 
sites. Frustration in this case refers to the impossibility 



of choosing the optimal phase difference for each bond. 
In a cold atom optical lattice system, frustration can be 
introduced by altering the phase factors for atom hop- 
ping between optical potential minima more explicitly, 
for example by following procedures similar to those pro- 
posed recently by Jaksch and Zoller , Mueller 0| , and 
S0rensen et al. The laser configurations suggested 

in these papers also enable spatially periodic modulation 
of the magnitude of boson hopping amplitudes, a feature 
that is important to the proposal outlined below. 

In a FF square-lattice JJA the sum of the optimal 
phase differences for individual bonds around every pla- 
quette is tt, fully incompatible with the integer multiple 
of 27r phase winding constraint imposed by the single- 
valued condensate wavefunction. For square lattice JJA's 
full frustration can be introduced by applying an external 
magnetic field that generates one half of a superconduct- 
ing flux quantum through each plaquette of the array. In 
the Landau gauge the frustration is imposed by changing 
the sign of every second vertical hopping parameter. For 
a FF-JJA, the Gross-Pitaevskii mean-field equation of 
the corresponding boson Hubbard model has two distinct 
degenerate solutions, illustrated schematically in Fig. ^ 
which break the discrete translational symmetry of the 
lattice, and for each solution a free overall phase factor in 
the condensate wavefunction which breaks gauge symme- 
try. The surprising property of FF square lattice JJA's, 
and by extension of FF square lattice cold atoms, is that 
the Ising order and the quasi-long-range phase order ap- 
pear to vanish nearly simultaneously and continuously at 
a common critical temperature. When quantum fluctu- 
ations are included, similar phase changes are expected 
to occur at zero temperature as the on-site interaction 
strength is increased. If these orders do in fact disap- 
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FIG. 1; Ground-state degenerate solutions for a classical FF- 
JJA. Double vertical lines stand for modulated "antiferro- 
magnetic" bonds (—a Ej) while single vertical and horizon- 
tal lines stand for unmodulated "ferromagnetic bonds" {E,j). 
The configuration shown corresponds to a = 0.5 for which 
Oa = ~0b = 7r/12 and Od = -Oc = 7r/4. The distinct 
configuration with equal energy is obtained by {9 a, do) 
— {9a,0d) or equivalently by vertical translation by one lat- 
tice constant. 

pear simultaneously, the phase transition would have to 
be in a new universality class and could not have a nat- 
ural description in terms of the condensate wavefunction 
order-parameter, a situation reminiscent of the decon- 
fined quantum critical behavior discussed recently [l^ 
by Senthil et al. 

The compound phase change in a frustrated JJA is 
closely related to the phase changes that occur in the vor- 
tex lattices of the mixed state of superconductors, and in 
rotating ^He and cold atom systems UU- The vortex 
lattice ground state has broken translational symmetry 
instigated by frustrating order-parameter-phase depen- 
dent terms in the Hamiltonian. The key difference be- 
tween vortex lattices and frustrated JJA's is that the 
broken translational symmetry is discrete rather than 
continuous in the latter case. Thermal fluctuations of 
a vortex lattice imply [isf that quasi-long-range phase 
order cannot exist at any finite temperature in 2D sys- 
tems. For superconductors it has been argued that 
given the absence of phase coherence, broken transla- 
tional symmetry will not occur either. For the FF-JJA 
case, the opposite conclusion has been reached in a care- 
ful Monte Carlo study by Olsson he finds that that 
vortex position fluctuations suppress the phase stiffness 
and instigate a BKT transition as the Ising phase tran- 
sition temperature is approached from below. If correct, 
this conclusion would have to be altered when frustration 



is weakened, as described below, and the Ising transition 
temperature is driven to zero. In this Letter we point 
out that these subtle phase changes can be studied by 
measuring the momentum distribution function (MDF) 
of a FF cold atom cloud, and report on theoretical esti- 
mates for the MDF based on a self-consistent harmonic 
approximation (SCHA) fiaj . 

We assume that atom hopping between sites on the 
optical lattice is weak enough to justify a single-band 
Wannier basis [if with Wannier function ii;(x). The lat- 
tice Hamiltonian we study is 

^/ = ^ E "x, - E cos (-^x. - k,+s) (1) 

where = d{n,m) with 71,111 £ [—J\f,J\f) is on a 2D 
square lattice with lattice constant d, S is the vector con- 
necting a lattice site to its neighbours, and the Joseph- 
son energy or atom hopping energies E:^, g are identi- 
cal (equal to Ej) on all bonds except the vertical bonds 
on every second column. These modulated frustrating 
bonds have the value —aEj with a > 0. In Eq. (P) 
the phase operator has been introduced by approx- 
imating the atom annihilation operator on site X; by 
— exp (i0xi ) I allowed when the mean occupa- 
tion n on each lattice site is large. The density fixi 
and phase (jj^i operators are canonically conjugate on 
each site. The negative hopping parameters introduce 
frustration, which can be energetically weakened [l^ by 
choosing a < 1. 

When quantum fluctuations are neglected, the T = 
condensate phase pattern is determined by minimizing 
the classical energy with respect to the phase difference 
X across positive Ej links; the single-valued condition 
requires that the magnitude of the phase difference across 
negative Ej links x' = ~3x, implying 's'l that sin(x) = 
Q;sin(3x) and hence that 

X = ± arcsin ( ^/[{3a- l)/a]/2) (2) 

for a > 1/3, while x = for a < 1/3. For a < 1/3, the 
energy penalty of frustration is paid completely on the 
negative Ej link and the classical ground ground state 
condensate phase is spatially constant. As a increases 
beyond this value, the energy penalty of frustration is 
increasingly shifted to the positive Ej links. The ground 
state configuration in this regime is doubly degenerate 
with currents circulating in opposite directions around 
alternating plaquettes, as illustrated in Fig. ^ Thermal 
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and quantum fluctuations will degrade both Ising and 
phase coherence orders. 

Phase coherence of cold atoms in an optical lattice 
can be directly detected by observing a multiple matter- 
wave interference pattern after ballistic expansion with 
all trapping potentials switched off. As time evolves, 
phase-coherent matter waves that are emitted from each 
lattice site overlap and interfere with each other. Narrow 
peaks appear in the MDF due a combination of lattice 
periodicity and long-range phase coherence 0, 0| . 
The vortex superlattice of the a > 1/3 mean- field state 
results in the appearance of additional peaks in the MDF; 
n/(k) = 3fie('I'^^(k)^'(k))/A where A is the system area, 
and \E'(k) is the 2D Fourier transform of the field opera- 
tor, 5'(x) = J2x "^(^ ^ ^i)bxi- It follows that 
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where we have defined a Wannier function form factor 
u'(k) = J cfxe^^^'^w{x) and the phase-phase correlator 
C(xi,Xj) = (exp [j(0xi — '/'xj)])- In the broken transla- 
tion symmetry state n f (k) is non-zero at superlattice re- 
ciprocal lattice vectors Gn,m = Tr{n,m)/d; for the classi- 
cal (i.e. U = 0) ground state at zero temperature we find 
that nf{G) = {N^ /A)n\w{G)\^ So{G) where N, = 4JV^ 
is the total number of lattice sites, and the superlattice 
structure factors are 



by that of an effective harmonic model defined by mean 
condensate phases on each site and harmonic coupling 
constants K on each nearest neighbour link. Minimizing 
the variational free-energy with respect to mean phases 
enforces average current conservation at each node of the 
lattice. Minimization with respect to the harmonic cou- 
pling constants sets them equal to the self-consistently 
determined mean curvature of the Josephson interaction. 
The phase changes across the vertical and horizontal pos- 
itive Ej links. Oh and 9y, are unequal in this approxi- 
mation, as are the harmonic coupling constants Kh and 
Ky and (of course) the coupling constant on frustrated 
links Ka. For [/ -> and T ^ 0, the dh = 9^ ^ x, 
Kh = Ky Ej cosx and — > ~aEj cos (3x)- 

The SCHA phase correlation function C(xi,Xj) = 
C^p Cq "^(Xi, Xj) is the product of a long-range factor 
C(xi,Xj) (3) C^p, dependent only on position within the 2x 2 broken- 
symmetry unit cell, and a Gaussian factor CQ''^(Xj,Xj) 
which captures the power-law decay of phase correlations 
in 2D superfiuids (here X^ is a lattice vector of the large 
unit cell so that sites are labelled by /i and i). We find 



that C^p is given by 
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(5) 



5*0(00,0) 
'S'o(Gi_o) 
'S'o(Go,i) 
50(01 1) 



[cos (x)cos (x/2)]^ 
[sin (x) sin (x/2)]^ 
[sin (x) cos (x/2)]2 
[cos (x) sin(x/2)]2 



and that 
C^'''(Xj,Xj) = cxp 



iV2 



EE 



(4) 



[l + 2iVBE(ek,./fcBr)] 



(6) 



with 5o(G„+2A;,m+2fc) = S'o(G„,„i) for any integers n, 
m, and k. Phase coherence in a lattice leads to con- 
densation peaks in n/(k) at all reciprocal lattice vec- 
tors G2„,2m- Coherence and Ising broken translational 
symmetry leads to additional peaks (satellites) with the 
characteristic pattern of structure factors summarized by 
Eqs. Q at the 2x2 superlattice reciprocal lattice vec- 
tors. MDF measurements therefore probe both types of 
order. 

These results will be altered by both quantum and 
thermal fluctuations. At low temperature (fc^T <C 
Ej) and well inside the superfluid regime {U <C Ej), 
the phase correlation functions are given reliably by a 
SCHA fli in which the density matrix is approximated 



In Eq. ® the sum is over the four Bogoliubov eigenmodes 
of the harmonic Josephson term at each wavevector in the 
2x2 super cell's Brillouin zone. Because we have cho- 
sen strictly on-site interactions, the quantum harmonic 
problem can be solved by first diagonalizing the Joseph- 
son interaction term, as in the classical case, and then 
performing independent Bogoliubov transformations on 
each mode. The contribution of a given Bogoliubov mode 
to the mean square phase difference between sites (/i, i) 
and (ly, ?' ) in Eq. © is therefore characterized by the 
quantity 

.Fi^;(X,-X,) = \v;ik)\' + \v:{k)f-2^e{y^{k)r 
X <(k)e*-(^'-^^+'^^'') } (7) 
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FIG. 2: The structure factor 5(k) for FF cold bosons in a 2D 
array with a — 0.5 as a function of the continuous variable 
kd e [-27r,27r] x [-2tv,2tt]. Here T = 0.242£;j//cs and 
)7 = O.lSj. 



where h^^, is the site separation for i — N-be{x) is 



a Bose-Einstein thermal factor, — t^'^k.o-, Ak,<j and 
w^(k) being the eigenvalues and the /i-th component of 
the eigenvectors of the harmonic Josephson interaction. 

We have evaluated S{k) = nf{k)A/{nN^\w{k)\'^) in 
the presence of both quantum and thermal fluctuations 
by summing over a finite lattice with Ng = 1296 sites 
in Eq. ||2Jl and applying periodic boundary conditions to 
make the wavevectors in Eq. © discrete. A typical re- 
sult is reported in Fig. [21 The presence of non-zero Ising 
satellites at k = Gi^o,Go,i and Gi^i is evident. These 
peaks are a sharp manifestation of the broken discrete 
translational symmetry and would be absent in an un- 
frustrated system. 

The evolution of S(k) with U at fixed T = 0.242 Ej/ks 
is illustrated in Fig.Owhere we plot 5(G) for Go,Oi Gi^o, 
Go,i and Gij. All four peaks are shghtly suppressed by 
quantum and thermal fluctuations with respect to the 
[/ = T = values in Eq. Q . At the critical value C/fg « 
0.14 the Ising satellites disappear while the condensation 
peak survives (the first order character of this transition 
is an artifact of the SCHA). The superlattice peaks may 
be regarded as Ising order parameters ^ S = sin (Oh) (see 
Eq. I0J). At the Ising point 9h 0, causing S{Go,o) to 
increase with increasing U, before resuming its decline. 

In summary, we have shown that FF cold atoms can 
offer a unique opportunity for experimental study of a 
system in which there is competition between critical 
phenomena associated with Z2 and gauge U(l) broken 
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FIG. 3: Condensation and Ising peaks of the structure factor 
S{G) as a function of U/Ej. The value of ^(G) for Gi,o, Go,i 
and Gi,i has been multiplied by a factor of 10 for clarity. The 
vertical dashed line indicates the value of Uf<^. 
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